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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

LASER  INDUCED  INTERMODULATED  FLAME  FLUORESCENCE: 
APPLICATION  TO  ANALYTICAL  ATOMIC  FLUORESCENCE 

By 

Lucas  P.  Hart,  III 
April  1984 

Chairman:  James  D.  Winefordner 
Major  Department:  Chemistry 

Many  transitions  used  in  atomic  flame  fluorescence  can 
be  saturated  with  modern  high  power  lasers.  The  fluores- 
cence signal  reaches  a constant  value  as  the  laser  irradi- 
ance  is  increased.  if  a cw  laser  is  modulated,  saturation 
of  the  transition  will  cause  nonlinear  response  to  the 
time-varying  excitation  waveform  and  generate  harmonics  in 
the  fluorescence  signal.  If  two  laser  beams  are  modulated 
at  different  frequencies  and  the  two  beams  are  combined,  sa- 
turation of  the  transition  will  produce  signals  at  the  in- 
termodulation product  frequencies.  Under  saturation  condi- 
tions, the  signals  at  these  frequencies  may  be  a significant 
fraction  of  the  signal  at  the  fundamental  frequency. 

Scatter  ,e.g.,  Rayleigh  or  Mie  scatter,  is  a linear 
process  and  does  not  produce  such  additional  frequencies. 


IV 


Measurement  at  the  intermodulation  product  frequency,  under 
saturation  conditions,  may  significantly  reduce  the  effect 
of  scatter  upon  the  fluorescence  signal  measured  at  a reso- 
nance line.  Additionally,  the  laser  beam  cross-sectional 

-3  2 

area  can  be  made  quite  small  (10  cm  ) so  that  the  intermo- 
dulation resulting  from  two  crossed,  modulated  beams  may  be 
a powerful  tool  for  spatial  diagnostics. 

The  intermodulation  signal  for  Na  in  a premixed,  lami- 
nar H2~02-Ar  flame  has  been  studied  using  a cw  dye  laser. 
Strong  intermodulation  signals  were  observed.  The  interfer- 
ence from  scatter  of  the  laser  beam  was  greatly  reduced  by 
measuring  fluorescence  at  one  of  the  intermodulation  product 
frequencies. 
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CHAPTER  ONE 
INTRODUCTION 


Laser  induced  fluorescence  spectroscopy  has  become  a 

valuable  analytical  tool  in  a wide  variety  of  fields,  now 

that  lasers  are  readily  available  with  wavelengths  in  the 

UV-visible  region  (1,2, 3, 4).  It  has  the  advantages  of  good 

sensitivity  and  spatial  resolution.  Lasers  may  be  focused 

to  beam  diameters  of  only  tens  of  micrometers,  giving  the 

potential  for  high  spatial  resolution  measurements,  or  the 

use  of  micro-samples.  Laser  induced  fluorescence  has  found 

use  in  ultra-trace  analysis,  where  the  goal  may  be 

single-atom  detection  (5).  Concentrations  as  low  as  10^ 
3 

atoms/cm  have  been  measured  for  sodium  in  atmospheric  pres- 
sure flames  (6).  Various  types  of  laser  induced  fluores- 
cence have  found  application  in  combustion  diagnostics 
(2,3,7),  where  techniques  originally  developed  using  atomic 
fluorescence  of  sodium  have  been  extended  to  the  measurement 
of  concentrations  of  various  radicals  throughout  the  combus- 
tion process. 

One  type  of  fluorescence  spectroscopy,  which  only  be- 
came feasible  with  the  availability  of  high  power  lasers,  is 
saturation  spectroscopy.  in  fluorescence,  the  signal  is 
proportional  to  the  excited  state  population.  The  excited 
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state  may  deactivate  either  r adiationally  or  by  collisional 
deactivation.  Thus  the  collisional  environment  may  affect 
the  fluorescence  signal.  However,  if  the  excitation  beam 
intensity  is  high  enough,  the  population  of  the  upper  level 
reaches  a maximum  as  the  rate  of  stimulated  emission  exceeds 
the  rate  of  spontaneous  emission  or  collisional  deactiva- 
tion. The  saturation,  as  this  condition  is  called,  may  be 
used  to  advantage  in  different  ways.  The  fluorescence  sig- 
nal is  determined  by  the  upper-level  population,  which  is 
now  independent  of  the  collisional  environment,  so  satura- 
tion may  be  of  value  in  diagnostic  studies,  where  the  colli- 
sional environment  may  change.  For  this  reason,  saturation 
spectroscopy  has  found  extensive  application  in  combustion 
diagnostics,  where  the  interest  is  in  radical  concentrations 
throughout  the  combustion  zone.  If  the  two  levels  (consid- 
ering a two-level  system)  have  equal  statistical  weights, 
the  upper  level  will  have  half  the  total  population  with 
full  saturation.  With  saturation,  the  fluorescence  signal 
for  a given  analyte  density  will  also  be  a maximum.  That 
is,  since  the  maximum  number  of  atoms  are  in  the  upper 
state,  the  fluorescence  transition  rate  will  also  be  maximum 
as  the  fluorescence  transition  cate  is  the  upper  level  den- 
sity times  the  spontaneous  transition  probability.  In  an- 
alytical atomic  spectroscopy,  saturation  spectroscopy  can  be 
applied  to  ultra-trace  analysis  and  in  the  quest  for 
single-atom  detection. 
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One  practical  limit  to  the  sensitivity  of  laser  induced 
fluorescence  spectroscopy  is  scatter  of  the  laser  radiation 
(8),  whether  it  be  Rayleigh  scatter  from  the  flame  gases, 
Mie  scatter  from  particulates,  or  scatter  from  the  sample 
cell.  Several  techniques  have  been  used  to  avoid  these  ef- 
fects. One  of  the  most  effective  is  to  use  non-resonance 
fluorescence  (1),  that  is,  to  look  at  a different  transition 
from  the  excited  state  to  a lower  state.  This  may  not  be 
practical,  due  to  a lack  of  other  transitions  or  a low  tran- 
sition probability.  Advantage  may  be  taken  of  the  polariza- 
tion of  Rayleigh  and  Mie  scatter,  and  polarizers  used  to 
reduce  the  scatter-to-s ignal  ratio  (8,9).  Two  light  sources 
may  be  used,  as  in  conventional  background  correction  tech- 
niques, but  with  a large  scatter  signal  and  small  fluores- 
cence signal,  statistical  uncertainties  may  be  large  as  the 
result  of  subtracting  two  large,  uncorrelated  numbers. 
Wavelength  modulation  of  a cw  dye  laser  has  been  used  to 
correct  for  scatter  resulting  in  an  order  of  magnitude  im- 
provement over  amplitude  modulation  (10) . 

Another  approach,  proposed  by  Alkemade  (5,11),  is  a 
variant  of  the  intermodulation  technique  originally  applied 
to  high-resolution  spectroscopy  by  Sorem  and  Schawlow  (12)  . 
In  this  approach,  the  laser  beam  is  split  into  two  beams  of 
equal  intensity,  each  is  amplitude-modulated  at  a different 
frequency,  and  then  the  beams  are  passed  through  the  atom 
cell  in  opposite  directions,  intersecting  therein. 
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Different  configurations  are  shown  in  Figure  1-1.  The 
scatter  will  be  proportional  to  the  sura  of  the  intensities 
of  the  two  beams.  Under  conditions  of  full  saturation,  the 
peak  fluorescence  signal  will  be  the  same,  independent  of 
whether  one  or  two  beams  are  present.  In  this  case, 
side-bands  will  be  generated  at  frequencies  equal  to  the  sum 
and  difference  of  the  two  modulation  frequencies.  This  ap- 
proach is  similar  to  one  evaluated  by  Freuholz  and  Gelbwachs 
(13),  in  which  the  laser  beam  is  sinusoidally  modulated; 
the  fluorescence  signal  also  appears  at  higher  harmonics  in 
the  case  of  saturation. 

In  this  work,  the  application  of  intermodulated  flu- 
orescence to  atomic  flame  spectroscopy  was  considered.  The 
system  used  for  experimental  characterization  was  a H2-02_Ar 
flame,  and  a cw  dye-laser  to  saturate  sodium  atomized  in  the 
flame.  The  behavior  of  sodium  atoms  in  flames  has  been  ex- 
tensively investigated  and  documented  (14) . Sodium  atoms 
are  readily  saturated  by  cw  dye  lasers.  The  cw  dye  laser  is 
appropriate  for  investigation  of  the  intermodulated  fluores- 
cence, as  it  is  stable  and  easily  modulated,  which  cannot  be 
said  of  many  higher  powered  pulsed  lasers. 

The  characteristics  of  various  optical  configurations, 
e.g.,  single-beam  modulation,  double-beam  modulation,  tight- 
ly-focused laser  beams,  or  beams  defocused  at  the  point  of 


Figure  1-1.  Optical  configuration  for  laser  induced  flame 
fluorescence  saturation  spectroscopy. 

a)  single-beam  configuration  for  harmonic  saturation 

b)  double-beam  configuration  with  nearly  co-linear 
beams  for  maximum  intermodulation  signal 

c)  double-beam  configuration  with  crossed  beams  for 
high  spatial-resolution  intermodulation  signal 


6 


flame 


o 

detector 


7 


observation,  were  studied  to  determine  the  range  of  applica- 
tion of  the  intermodul a ted  technique.  Real-time  measure- 
ments using  a phase-sensitive  amplifier  (lock-in  amplifier) 
permitted  the  optimization  of  various  parameters  and  were 
used  to  demonstrate  the  application  of  intermodulated  flu- 
orescence to  analytical  measurements.  The  fluorescence  sig- 
nals were  also  digitized,  and  power  spectra  calculated  using 
a fast  Fourier  transform  ( F FT)  to  simultaneously  determine 
the  magnitudes  of  the  signals  at  the  various  intermodulation 
products  frequencies,  in  order  to  characterize  the  behavior 
of  the  intermodulated  fluorescence  signal. 

A qualitative  description  of  the  behavior  of  the  inter- 
modulated fluorescence  was  derived  using  a basic  model  of  a 
two-level  atomic  system,  the  rate  equation  approach  to  solv- 
ing the  steady-state  system  response,  and  Fourier  series  ex- 
pansion of  the  time-varying  modulation  waveforms. 

The  significant  reduction  of  scatter  resulting  from  the 
use  of  an  intermodulation  product  frequency  for  the  fluores- 
cence measurment  was  shown.  An  alternative  application  of 
intermodulated  flame  fluorescence,  for  spatial  resolution. 


was  demonstrated. 


CHAPTER  TWO 

THEORETICAL  CONSIDERATIONS 


The  theory  of  saturation  spectroscopy  has  been  well  do- 
cumented in  detail  for  a wide  range  of  applications,  includ- 
ing high  resolution  atomic  and  molecular  spectroscopy  (15)  , 
optical  pumping  (16),  and  analytical  spectroscopy  (1).  The 
intermodul a ted  absorption  signal  has  been  treated,  using 
several  different  formalisms.  Freuholz  and  Gelbwachs  have 
considered  the  case  of  single  beam  saturating  atomic  flu- 
orescence in  a flame  (13).  For  this  work,  a simplified, 
i.e.,  two-level,  system  will  be  considered  to  show  the  ori- 
gin of  the  intermodulated  fluorescence  signal,  first  for  a 
single-beam  system  and  then  for  a crossed-beam  system.  The 
two-level  system  is  adequate  to  show  the  origin  of  the  in- 
termodulated component  and  show  the  relationship  to  the 
laser  excitation.  A quantitative  analysis  would  have  to 
take  into  account  factors  such  as  the  Gaussian  intensity 
distribution  of  the  laser,  laser  induced  chemistry  and  ioni- 
zation, multiple  electronic  levels,  etc.  The  rate  equation 
approach  used  is  appropriate  for  atmospheric  pressure  flames 
(17)  . 

The  case  considered  is  an  atom  with  two  energy  levels 
that  is  a trace  constituent  in  an  atmospheric  pressure  sys- 
tem where  the  collision  rates  are  high  enough  for  coherent 
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effects  to  be  disregarded.  The  excitation  source  is  consi- 
dered to  have  a spectral  profile  much  greater  than  the  ab- 
sorption profile  of  the  absorbing  atoms,  and  thus  has  a uni- 
form spectral  density  and  interacts  with  all  the  absorbing 
atoms.  The  source  is  also  assumed  to  have  a uniform  spatial 
distribution.  The  two  levels  of  the  atomic  system  are  as- 
sumed to  have  equal  statistical  weights.  The  source  is 
tuned  to  the  atomic  transition,  with  a spectral  bandwidth 
Sv  (Hz)  , a power  $ (W)  , and  a cross-section  S (cm2)  . The 
population  of  the  levels  (cm  ) is  indicated  as  n^  and  n^ 

and  the  total  atom  population  ,n.  + n„ , as  n . 

1 2 t 

Using  a rate  equation  approach,  and  ignoring  collision- 
al  excitation,  as  shown  in  Fig.  2.1a,  the  population  of  the 
two  levels  is  given  by 

nl[B12Ev/cl  = VB2lVC  + A21  + k21l  ' (2a) 

where : 

B = (c3/8TThv3T  ) : Einstein  coefficient  for  stimulated 

absorptionpand  emission,  J cur  Hz; 

A = Einstein  coefficient  for  spontaneous  emission,  s_1; 

k = collisional  rate  constant,  s”  ; 

c = light  velocity  in  vacuum,  cm  s-1; 

h = Planck  constant,  J s; 

v = resonance  frequency  of  the  transition,  Hz; 

t = spontaneous  (radiative)  lifetime  of  the  upper  level,  s 

SP  ( = 1/A) ; 

t = lifetime  of  upper  level  (including  collisions)  s 
(=  1/A+k)  ; 
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Ev>  = spectral, irradiance  of  the  laser  beam, 

6vS  J s 1 cm  £ Hz  ; 

The  saturation  parameter  or  saturation  irradiance,  ES  , is 
defined  as  that  irradiance  which  would  give  an  upper-level 
population  half  that  of  an  infinite  irradiance  (1).  Then, 

s 3 

Ev  = 8 hy  tso  = saturation  parameter  for  the  atomic 
2C2  x P transition,  J s cni  Hz-1. 

Substituting  both  the  saturation  parameter,  , and  nT-n2 
for  n.^  in  Eqn.  2.1  gives: 


n2  1 


K 


(2.2) 


nT  2 <+ 


This  is  the  steady-state  approximation;  the  relative  popu- 
lations of  the  two  levels  are  shown  in  Fig.  2.1b.  if  the 
atomic  relaxation  time  is  fast  compared  to  the  modulation 
rate,  the  time-varying  excited-state  population  is  given  by 


n„(t)  1 e (t) 

^ V 

nT  2 ES+  Ev  (t) 


(2.3a) 


or 


n2(t)  _ 1 E°f(t) 
nT  2 ES+  E°f(t) 

where  E^(t)  is  the  time  varying  spectral  irradianc 
peak  spectral  irradiance,  and  f(t)  the  function  d 
the  time  variation.  The  fluorescence  radiance  is 
Bp  = n2hv  U/4tt)A21 


(2.3b) 

o 

e , Ev  the 
escr ib ing 

(2.4) 


where  £(cm)  is  the  depth  of  the  fluorescence  volume.  The 


Figure  2-1.  Two-level  system 

a)  energy-level  diagram 

B = Einstein  coefficient  for  stimulated  absorption 
and  emission 

&2\  = Einstein  coefficient  for  spontaneous 
emission 

k21  = collisional  quenching  rate  constant 

b)  Saturation  curve  for  two-level  system 
n2,  n = population  density  of  levels  two 

and  one 

E = source  spectral  irradiance 

S ... 

Ey  = saturation  irradiance 
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E/Es 
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time-varying  fluorescence  signal  is  thus  proportional  to 
n2(t).  The  effect  of  varying  source  irradiance  will  be 
examined  for  single-beam  excitation,  with  arbitrary  wave- 
form, and  crossed  or  co— linear  beams  with  square— wave  or 
sinusoidal  modulation. 


Single-Beam  Modulation 

Considering  the  amplitude-modulated  excitation  waveform 
as  a 100%  modulated  positive  periodic  function,  it  may  be 
expressed  as 


Ev  (t)  = fv  [1+f  ( t)  ] 


(2.5) 


where  Ev  is  the  peak  irradiance  and  f(t)  could  be  sin(cot) 
for  sinusoidal  modulation  for  example.  To  simplify  the  no- 
tation in  the  following  treatments,  let  y(t)  = 2(n2/n  ). 
Then  the  equation  (2.3b)  may  be  rewritten  as 


y (t)  = 


or 


y(t)  = 


Ev  [1+f  ( t)  ] 
2ES  + E° [ 1+f ( t ) ] 

V V 


Ev  tl+f(t)] 

2EV  + EV  + 


(2.6) 


(2.7) 


Since  2ES  + E°+  E°f(t)  is  equivalent  to 

V V V ^ 


[2EvS  + Ev°][l  + 


V 


2ES+  E° 

V V 


f (t)  ] 


(2.8) 


equation  2.7  may  be  written  as 
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y(t) 


E°[l+f (t) 1 
2Ev  + Ev 


-] 


U+ 

2Ev+  Ev 


(2.9) 


Let  a-  E^/2E^  + E^  and  z - af  (t)  . Since  |z[<l,  a power  series 
expansion  may  be  used  (18),  i.e., 

(1+z)  1 = l-  z + z2-z3  + (2.10) 

and 

Ev[1  + f(t)]  2 3 

y(t)  = ( i 5 ) (l-z+z  -z  + ....)  (2.11) 

2E  + E 

V V 

The  effect  of  saturation  upon  the  response  of  y(t)  to  the 
modulation  waveform  may  be  examined  by  substituting  for 


f (t)  . 


Sinusoidal  Modulation 

In  this  case,  f(t)  may  be  represented  as  either  sin(wt) 
or  cos(tot).  Keeping  terms  with  frequency  up  to  a power  of 
three 

E° 

y(t)  = ( F — o){1  + d-ct)  [f  (t)-af2(t)  +a2f3  ( t)  ] } (2.12a) 

2 E + E 

v v 

g ives 

o 

g 

Y(t)  = — sT — o{1  + (1-a)  tcoswt  - -a(l+cos2ait)  (2.12b) 

2E  + E . _ 2 

v v 12, 

+ ja  (cosu)t+cos3cot)  ] } 


It  can  be  seen  that  the  effect  of  increasing  a is  to  de- 
crease the  fraction  of  the  total  signal  with  a frequency 
component  the  same  as  the  modulation  frequency,  and  to  in- 
crease the  harmonic  content. 
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Square-Wave  Modulation 

A square  wave  may  be  represented  as  a Fourier  series 

(19) 

2 1 1 

f(t)  = — (cosoit  - ycos  3a)t  + jcos  5wt  ) (2.13) 

From  (2.12a),  it  can  be  seen  that  the  frequency  components 
of  the  original  excitation  waveform  are  attenuated.  To  de- 
termine the  resultant  frequency  components,  the  second  power 

2 

term  may  be  examined.  The  (l-a)af  (t)  may  be  expanded, 
keeping  terms  up  to  the  5th  harmonic 

f 2 ( t)  = [2  (cosoot  - ycos3oot  + ycos5wt)]2  (2.14a) 

4 2 1 1 

= -2  (cos  cut  - ycos3iotcos  mt  + — gos5mtcoswt 

- ycoswtcos3 wt  + ycos3  wtcos3  tot  - yycos5wtcos3 cut 

+ ycosu}tcos5o)t  - yycos3 wtcos5 mt  + yycos5  wtcos5mt) 

(2.14b) 

Since  cos(A)cos(B)  = 1/2  ( cos  (A+B) +co  s ( A- B)  ),  this  may  be 
expressed  as 

2 2 1 1 
f (t)  = — [l+cos2u)t  -rr(cos2mt+cos4cot)  +—  (cos4^ t+cos6^ t) 

7T  J 5 

-y  (cos2wt+cos4(jot)  +y(l+cos6mt)  -yy  (cos2mt+cos8^t) 

+y  ( CO  S4(0  t+CO  S6t0  t)  — yy  ( CO  S2u)  t+CO  S8  U)t ) +yy  ( 1 + COSlOOJt) 

(2.14c) 

Continuing  the  expansion  of  this  series  gives  coefficients 
°f  If  -1/3,  -1/15,  -1/35,  -1/63,  -1/99  or  a second  harmonic 
term 

f2(t)  = ( 1/tt2)  cos2a)t 
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Since  2.14c  contains  only  even  terras,  it  can  be  seen  that 
the  the  third-power  term  will  have  no  second  harmonic  terra, 
while  the  second-harmonic  terras  from  the  fourth-power  terra 
have  a coefficient  of  a2 /4  and  will  be  smaller.  The  ratio 
of  the  second  harmonic  to  the  fundamental  is  then  approxi- 
mately a/2ir. 


Double-Beam  Modulation 

The  effect  of  modulating  two  intersecting  beams  at  dif- 
ferent frequencies  may  be  examined  using  the  same  approach 
as  for  the  single-beam  case,  but  first,  a heuristic  approach 
to  the  square-wave  modulation  more  clearly  illustrates  the 
effects,  and  gives  an  estimate  of  the  coefficients  for  the 
intermodulated  terms. 

Square-Wave  Modulation 

In  Figure  2-2,  the  first  two  waveforms  are  those  of  the 
excitation  waveforms.  If  there  is  linear  combination,  or 
superposition,  of  the  two  waves,  the  amplitude  will  be  the 
sum  of  the  two,  giving  a result  as  in  Fig.  2.2c.  In  the 
case  of  intermodulation,  or  nonlinear  response,  the  ampli- 
tude of  the  resulting  waveform  may  be  less  than  the  sum  of 
the  two  waves,  as  in  Fig  2 . 2d . The  analysis  of  the  intermo- 
dulation process  is  considerably  simplified  by  considering 
the  resultant  waveform  as  being  the  superposition  of  the 
train  of  square  waves  in  (a)  and  modulated  train  of  waves 
with  frequency  of  (b) . The  modulated  train  may  be  expressed 


Figure  2-2.  Waveforms  for  intermodulated  f luorescence, 

a)  square-wave  modulation  of  beam  1,  amplitude  2Yq 

b)  square-wave  modulation  of  beam  2,  amplitude  y 

c)  linear  combination  of  a and  b waveforms 

d)  saturation,  non-linear  combination  of  a and  b 

e)  modulation  waveform,  when  multiplied  by  waveform 
in  b and  added  to  waveform  in  a,  gives  d 
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as  the  product  of  the  train  in  (b)  , y(t),  and  a modulation 
function , M(  t)  , i .e . , 

1 + 2y  ( t)  = 1y(t)  + 2y  ( t)  M ( t)  (2.15) 

The  modulation  function,  M(t)  shown  in  (e),  has  the 

same  frequency  as  the  wavetrain  in  (a),  and  an  amplitude 

varying  between  1 and  y , where 

J m 


y 


m 


1+2  o 1 o 

y - y 


2 o 

y 


and  '''+2y0  is  the  amplitude  when  both  beams  are  combined. 
For  the  case  of  square-wave  modulation  being  considered,  the 
time-varying  functions  may  be  expressed  as 


y(t)  = | [l+f1  ( t)  ] 

9 2 0 

y (t)  = J [1  + f 2 ( t)  ] 


(2.16a) 

(2.16b) 


The  modulation  function  may  be  expressed  as 
1+2  o 1 o 1+2  o 1 o 

M ( t)  = t[(1+  Y " Y + (1-  Y " Y ) f,  (t)]  (2.16c) 

2 2o  2 o 1 

y y 


The  square-wave  function  may  be  represented  by  the  Fourier 

series  expansion  (eqn.  2.13).  Keeping  only  those  terms 

with  the  fundamental  frequency  in  each  term,  we  get  for  the 

intermodulation  term,  y.  (t) 

J lm 

_ 2.1  o^2o  1+2  o, 

Yim(t)  -TT“(  y + y - y ) cos  u tcos  co  t (2.17) 

The  height  of  the  square  wave,  with  a given  excitation, 
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ie  1 o 2 o o 2 o 

y y / or  ly°+  y°  raay  be  determined  from  the  basic 

relationship,  Eqn.  2.2.  After  algebraic  manipulation,  the 

amplitude  of  intermodulation  product  may  be  expressed  in 

terms  of  the  excitation  amplitudes  (y.  (t)  = C.  f(t)  ) as 

lm  lm 


1E° 

v 


2E° 

v 


We° 

v v , 


Clm  ’‘\°+E.S  + 2E°+ES  ' 1E°+2E°+E=) 


(2.18) 


v v 


V V V 


From  examination  of  this  equation,  it  is  apparent  that  at 

low  intensities,  i.e.,  ^E°  or  2e°  <<  E^,  the  inter  modulated 

term  goes  to  zero.  At  high  intensities,  1E°,2E°  >>  EvS,  C. 

2 

goes  to  1/tt  . For  a given  laser  power,  the  maximum  intermo- 
dulation signal  is  obtained  if  the  power  is  divided  equally 
between  the  two  beams,  in  which  case 


C,-m  = -A 


(E°)2 


im  it2‘(E°+ES)  (2E°+ES) 

V V V V 


(2.19) 


where  the  peak  intensity  E°  is  the  same  for  either  beam, 
o 1 o 2 o 

i.e.,  E^  = Ev  = Since  the  component  at  either  funda- 
mental,in  the  lack  of  the  other  beam  being  on,  is  given  by 
equations  2.13  and  2.16,  the  ratio  of  the  component  at  the 
intermodulation  frequency  to  that  at  the  fundamental  (one 
beam  blocked  so  there  is  no  intermodulation)  is 


Cim 


v 


2E°  + ES 

V V 


(2.20) 


It  should  be  noted  that  this  is  a first-order  approximation. 
Intermodulation  terms  as  a result  of  the  fundamental  and 
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third  harmonic  term  of  the  Fourier  series  expansion  of  the 
modulation  waveform  will  give  components  at  (3u^+  co2)  / 

( 3 a)  2)  / ( w^+3  m 2)  , and  ( 3 co  - oj^),  each  with  a coefficient 
one-sixth  that  of  the  intermodulation  term  resulting  from 
the  fundamentals  only.  This  heuristic  approach  gives  a sim- 
ple presentation  based  upon  square-wave  modulation.  For  an 
arbitrary  waveform,  an  analysis  similar  to  that  for  the  sin- 
gle-beam system  can  be  performed. 

Arbitrary  Waveform  Modulation 

For  an  arbitrary  waveform,  the  basic  equations  used  for 
the  single-beam  case  may  be  expanded  for  the  case  of  two 
crossed  or  nearly  co-linear  beams  modulated  at  different 
frequencies.  The  instantaneous  intensity  will  be  the  linear 
sum  of  the  intensities  of  each  of  the  beams  and  Eqn.  2.5 
becomes 

2E° 

E (t)  = — [1+  f,  (t)  ] + — — [1  + f (t)  ] (2.21) 

2 1 2 1 

These  equations  only  apply  to  100%  modulation  of  each  beam. 
A more  general  treatment  would  include  a possible  offset  in 
the  dc  term  as  was  done  for  the  modulation  waveform  in  the 
two-beam  square  treatment  (Eqn.  2.16c).  For  clarity  and 
ease  of  notation,  only  100%  modulation  will  be  discussed. 
Following  the  same  algebraic  steps  used  in  the  single-beam 
case,  the  time-varying  population  of  the  upper  level  may  be 
wr itten 
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y(t)  = 


1E°[i+f1(t)] 


+ 2E°[l+f2(t)] 


[2ES+1E°+2E°]  [1+- 

V V V 


s. l„o' 2„ofl ^ ^ 


2eu+-e-+-ev 


2es+1e°+^e° 

V V V 


(2.22) 

-f,(t)  ] 


This  is  similar  to  Eq.  2.9.  Using  the  series  expansion, 
where  now 


z = 


V 

rxr 


2ES+1E°+2E° 

V V V 


-f1  (t) 


2ES+ 

v 


V 

Wi 

v 


*f2(t) 


(2.23) 


we  obtain  a general  equation  for  the  two-beam  case,  modulat- 
ed at  different  frequencies 


y(t)  = ~l  o12  oiVtl+fiUH+Vtl+f  (t)]}[l-z+z2+...] 

2 E + E + E vi  v Z 

V V V 

(2.24) 

To  compare  this  to  the  previous  treatment  of  the  square-wave 
modulated  two-beam  system,  consider  the  fundamental  frequen- 
cy only  of  the  square  wave  as  before,  the  first  power  term 
of  the  power  series  expansion,  and  let  E0=1E0=2E°.  Equation 

V V V 

2.24  then  simplifies  to 
E°  - 

Y(t)  = ^[1  +^(cos(i)1t+cos^t)][l-z]  (2.25) 

E + E 

V V 

which  gives  an  intermodulation  coefficient  of 


Cim  = f ( 

7T 


V 


s o 
E + E 

V V 


(2.26) 
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This  result  is  quite  similar  to  that  obtained  in  Eq.  2.19. 
Thus  both  approaches  give  the  same  qualitative  results, 
showing  how  the  intermodulation  component  originates.  The 
first  treatment  of  square-wave  modulation  for  the  two-beam 
system  is  conceptually  simpler;  the  more  general  approach, 
based  on  an  arbitrary  waveform,  can  be  extended  to  compare 
the  effects  of  varying  degree  of  modulation,  different  wave- 
forms, and  real  conditions  which  have  been  avoided  in  the 
simplified  model  used  here,  e.g.,  the  effects  of  Gaussian 
spatial  intensity  distribution  could  be  included  and  the  ap- 
propriate equations  numerically  solved. 


CHAPTER  THREE 

EXPERIMENTAL  EQUIPMENT  AND  METHODS 

The  experimental  system  used  for  fluorescence  measure- 
ments is  schematically  illustrated  in  Figure  3-1.  It  may  be 
functionally  broken  down  into  the  optical  system,  for  excit- 
ing and  detecting  the  fluorescence,  the  atom  cell,  i.e.,  the 
flame,  and  the  electronic  detection  and  signal  processing 
equipment.  These  will  be  discussed  separately,  and  then  the 
experimental  procedures  described. 

Optical  System 

The  principal  optical  system  components  and  their 
manufacturers  are  given  in  Table  3-1.  A continuous-wave  dye 
laser,  pumped  by  an  argon-ion  laser,  was  used  as  the  excita- 
tion source.  The  theory  of  dye  lasers  and  their  application 
to  analytical  spectroscopy  has  been  extensively  discussed  in 
numerous  review  articles  and  books  (1,20,21,22);  only  the 
relevant  characteristics  will  be  discussed  here.  The  Ar-ion 
laser  has  a nominal  maximum  output  power  of  18  W when 
operated  in  the  all-lines  mode.  In  this  mode,  a broad-band 
mirror  is  used  as  the  rear  reflector  of  the  laser  cavity, 
and  lasing  occurs  simultaneously  at  the  wavelengths  of  sev- 
eral different  argon  ion  lines,  primarily  those  with 
wavelengths  of  514.5  nm  and  488.0  nm.  An  alternative, 
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Table  3-1.  Optical  Equipment  used  in  Intermodulated 
Fluorescence  Studies  and  its  Source  of  Manufacture 


Item 

Argon  ion  laser, 

Model  171-19 

CW  dye  laser.  Model  375 

Power  meter,  Model  210 

Beamsplitter,  cube 

Lenses,  spherical 
achromatic  doublets 

Lens  positioner,  5-axis, 
Model  LP-1 

Lens,  spherical  biconcex 

Polarizers,  Gian  Thompson 

Polarizer  rotators. 

Model  1641 

Chopper,  Model  7500 
Chopper,  Model  125 

Flame,  burner  head 

Nebulizer, mixing 
chamber.  Model  303 

Monochromator,  Model  H-10 


Source 

Spec tra-Phy sics 
Mountain  View,  Calif. 

Spectra-physics 

Coherent  Radiation 
Palo  Alto,  Calif. 

Melles  Griot 
Irvine,  Calif. 

Melles  Griot 


Newport  Corporation 
Fountain  Valley,  Calif. 

ESCO  Products,  Inc. 

Oak  Ridge,  N.J. 

Karl  Lambrecht  Corp. 
Chicago,  ill. 

Oriel  Corporation 
Stamford,  Conn. 

Rofin-Math  Associates 
Great  Neck,  N.Y. 

Princeton  Applied  Research 
Corp.,  Princeton,  N.J 

labor  a tor y-constr uc ted 

Perk  in-Elmer 
Norwalk,  Conn. 

American  ISA  Inc. 

Metuchen,  N.J. 


Photomultiplier,  Model  R777  Hamamatsu 

Middlesex,  N.J. 
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single-line  operation,  uses  a wavelength-selective  element, 
namely  a prism  in  front  of  the  rear  reflector.  Higher  out- 
put powers  are  obtained  for  a single  line,  than  in  the 
all-lines  mode,  but  the  total  power  is  less.  For  pumping 
Rhodamine  6G,  the  aye  used  here,  the  all-lines  mode  gives 
more  efficient  conversion  of  argon  ion  laser  pump  power  to 
dye  laser  output  power.  The  power  of  the  Ar-ion  laser  was 
monitored  with  the  power  meter  included  with  the  power  sup- 
ply. The  power  supply  has  two  possible  modes  of  operation: 
constant  current  or  constant  output,  in  which  a feedback 
circuit  is  used  to  vary  the  current  to  maintain  a constant 
output  power.  Adequate  stability  of  the  dye  laser  output 
was  maintained  using  the  constant  current  control  and  the 
constant  output  mode  was  not  used. 

The  continuous-wave  dye  laser  used  had  been  modified 
with  the  extended  cavity  option  (Spectr a-Physics  Model  341) 
for  alternative  use  as  a source  of  picosecond  duration  laser 
pulses.  The  output  mirror  for  that  option  is  not  optimum 
for  use  for  cw  excitation,  having  a different  focal  length 
and  transmission.  Therefore,  the  conversion  efficiency  was 
lower  than  would  normally  be  obtained  with  Rhodamine  6G. 
The  argon-ion  laser  was  operated  to  produce  9 W pump  power, 
all-lines,  which  produced  1.4  W at  the  589.6  nm  wavelength 
of  the  soaium  absorption  line,  using  approximately  3 
mM  Rhodamine  6G  (Exciton  Rhodamine  590)  in  reagent-grade 
ethylene  glycol.  Typically,  Rhodamine  6G  has  a conversion 
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efficiency  of  25%.  The  Lyot  or  birefringent  filter  used  as 
a wavelength  tuning  element  in  the  dye  laser  results  in  a 
linewidth  of  approximately  80  pm  (manufacturer's  litera- 
ture) , which  is  large  compared  to  the  typical  atomic  absorp- 
tion linewidths  in  flames  of  5 pm  (23),  so  that  this  laser 
may  be  considered  a quasi-continuum  source.  The  dye  laser 
uses  a free-flowing  dye  jet  resulting  from  pumping  the  dye 
solution  through  a small  nozzle  at  high  pressures.  The  re- 
commended operating  pressure  is  that  at  which  the  power  lev- 
els off  as  the  pressure  is  increased.  For  this  work,  a 
pressure  of  100  psig  was  used.  (Pressures  somewhat  higher 
than  that,  e.g.,  above  120  psig,  may  cause  the  magnetically 
coupled  pump  to  decouple,  resulting  in  loss  of  pressure  and 
possible  splattering  of  dye  over  the  optics,  or  the  nozzle 
to  come  free  of  its  holder,  ending  up  in  the  dye  reservoir, 
and  resultant  loss  of  lasing.) 

Optical  Components 

Low-loss  dielectric  mirrors  were  used  for  steering  the 
beam  from  the  cw  dye  laser  to  the  flame  system.  As  indicat- 
ed in  Figure  3-1,  a beamsplitter  was  used  for  measurements 
involving  two  beams.  The  beamsplitter  was  a cube  type 
beamsplitter,  composed  of  two  right  angle  prisms,  with  a 
metal-dielectric  film  deposited  on  the  interface  before  they 
are  cemented  together.  The  percentage  transmittance  and  re- 
flectance are  specified  as  45  and  55,  respectively; 
however , this  does  not  account  for  losses  as  a result  of  the 
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metal-dielectric  film.  The  transmittance  and  reflectance 
were  measured  to  be  37%  and  53%.  Single  beam  measurements 
were  made  with  the  beamsplitter  removed  and  only  the  beam 
path  on  the  right  of  the  diagram  (M2,  M3,  L2 , CH2)  used. 

For  all  dye  laser  output  power  measurements,  a Coherent 
absorption  type  power  meter,  with  a sensitivity  of  1 mV/W 
was  used,  in  conjunction  with  the  Coherent  readout,  or  a mi- 
crovoltmeter (Boonton  model  95A)  and  recorder  readout. 

Achromatic  doublet  lenses  with  a focal  length  of  300  mm 
were  used  to  focus  the  laser  into  the  flame.  Assuming  a 
Gaussian  spatial  distribution  and  the  measured  divergence  of 
2 mr ad  for  the  dye  laser,  the  cross-sectional  area  of  the 

waist  of  the  focused  beam  is  calculated  to  be  of  the  order 

-3  2 

of  10  cm  . Because  of  the  critical  positioning  required 
for  alignment  of  the  two  beams  with  respect  to  each  other, 
five-axis  lens  positioners,  with  fine  adjustment  screws  al- 
lowing better  than  1 ym  resolution,  were  used.  The  choppers 
were  placed  as  close  to  the  flame  as  space  would  permit  so 
that  the  wave  form  of  the  modulated  laser  beam  would  be  as 
close  to  a square  wave.  The  Rofin  chopper  (CHI)  was  operat- 
ed at  178  Hz,  and  had  a 9 mm  aperture  in  the  rotational 
direction.  The  PAR  chopper  (CH2)  was  operated  at  260  Hz  and 
had  a 15  mm  aperture  in  the  rotational  direction, 
blades  with  fan-blade  shaped  apertures.) 


(Both  had 
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Flame  System 

A nearly-stoichiometr  ic  H2~02'’Ar  fiame  was  supported  on 
a laboratory  constructed  Meeker  type  burner,  which  has  been 
used  extensively  in  this  laboratory  for  flame  diagnostic 
studies  (24).  The  brass  burner  head  has  96  holes,  in  a 8x12 
rectangular  grid,  32  for  the  inner  flame  in  which  the  sample 
is  atomized,  and  64  for  the  outer  flame  which  served  to  keep 
the  inner  flame  composition  and  temperature  uniform.  The 
inner  flame  cross  section  was  approximately  4 mm  by  8 mm. 
The  burner  head  was  mounted  on  a standard  premixing  chamber 
and  pneumatic  nebulizer  . The  gas  flow  rates  were  estimated 
to  be  4.5,  2.3,  and  9 L/min  of  hydrogen,  oxygen  and  argon 
respectively.  Another  2 L/min  of  argon  were  required  for 
the  nebulizer,  which  aspirated  2 mL/min.  All  measurements 
were  made  with  the  optical  axis  2 cm  above  the  burner  top. 
All  fluorescence  measurements  were  made  using  a sodium  con- 
centration of  1 yg/mL. 

For  a similar  flame,  van  Calcar,  et  al . (23)  calculat- 
ed and  measured  the  saturation  parameter  for  sodium;  they 
measured  a value  of  5.6x10  cm  2hz  ^ , in  terms  of  spec- 
tral volume  density,  ( p^s) , which  differed  from  the  calcu- 
lated value  by  less  than  the  experimental  error.  For  the 
laser  used  here,  with  a spectral  bandwidth  of  60  gHz,  that 
value  corresponds  to  a saturation  parameter,  in  spectral  ir- 
radiance,  of  about  100  W cm-2  (e®  = p®c6v  ).  Using  a value 
of  10  cm  for  the  cross-sectional  area  of  the  beam  waist, 
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a maximum  irradiance  of  1.2  kW  cm  can  be  calculated; 

thus,  the  transition  should  be  readily  saturated. 

Optical  Detection  System 

Fluorescence  was  collected  by  a 2 in  diameter,  2-1/2  in 
focal  length  lens  (L3)  and  focused  on  the  entrance  slit  of 
the  high  luminosity  f/3.5  monochromator.  An  aperture  was 
placed  in  front  of  the  lens  to  reduce  stray  light.  The 
monochromator  was  oriented  vertically,  with  the  slits  paral- 
lel to  the  laser  beam,  to  increase  the  volume  of  the  inter- 
secting laser  beams  that  was  viewed.  The  mimimum  slit 
height  of  2 mm  was  used.  The  0.5  mm  slit  width  used  gave  a 
spectral  bandwidth  of  4 nm , wide  enough  so  both,  the  589  .0 
and  589.6  nm  sodium  lines  were  detected.  The  photomulti- 
plier was  operated  at  -600  V to  keep  the  output  current 
level  below  10  yA  and  avoid  nonlinear  response  of  the  pho- 
tomultiplier tube. 

Electronic  Measurement  System 
Two  different  systems  were  used  for  measurements;  the 
components  are  given  in  Table  3-2.  A lock-in  based  system 
was  used  for  real-time  measurements.  The  results  of  making 
adjustments  in  the  optical  path,  tuning  the  laser,  etc. 
were  followed  using  this  system  as  were  quantitative  meas- 
urements after  the  system  was  adjusted.  The  computer-based 
digital  data  acquisition  system  was  used  to  measure  power 
spectra,  so  that  the  amplitudes  of  signals  at  the  various 
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Table  3-2.  Electronic  Equipment  used  in  Intermodulated 
Fluorescence  Studies  and  Its  Source  of  Manufacture 


Item 

High  voltage  power  supply, 
Model  226 

Preamplifier,  Model  113 

Lock-in  amplifier, 

Dynatrac  391 

Lock-in  reference 

Recorder,  dual  pen, 

Series  5000 

Anti-aliasing  filter 

Minicomputer,  PDP  11/34A 

Analog/digital  converter, 
Model  LPSAD-12 

Video  graphics  display, 
GT600  option  for  DEC  VT100 

Printer,  Model  IDS-445 

Plotter,  HiPlot  DMP-4 
Voltmeter,  DC,  Model  95A 


Source 

Pacific  Precision  Measurements 
Concord,  Calif. 

Princeton  Applied  Research 
Corp.,  Princeton,  N.J. 

Ithaco  Inc. 

Ithaca,  N.Y. 

laboratory- cons true ted 
(Appendix  A) 

Houston  Instruments 
Austin,  Texas 

laboratory-constructed 
(Appendix  A) 

Digital  Equipment  Corporation 
Maynard,  Mass. 

Digital  Equipment  Corporation 


MATROX  Electronic  Systems  Ltd. 
Quebec,  Canada 

Integral  Data  Systems,  Inc. 
Milford,  n.H. 

Houston  Instruments 

Boonton  Electronics  Corp. 
Parsippany,  N.J. 
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harmonics  and  intermodulation  frequencies  could  be  deter- 
mined simultaneously.  The  analog  system  was  used  to  monitor 
the  signal  to  make  sure  there  was  no  unusual  change  in  sig- 
nal during  the  digital  data  acquisition  period.  The  voltage 
signal  measured  by  each  was  generated  by  passing  the  current 
from  the  photomultiplier  through  a lOkfl  load  resistor.  That 
voltage  was  amplified  by  the  preamplifier,  normally  with  a 
gain  between  10  and  200.  The  lower  cutoff  frequency  of  the 
preamplifier  was  set  at  0.03  Hz  to  avoid  DC  drift  problems, 
and  the  upper  cut-off  frequency  set  at  3 kHz  to  attenuate 
noise  and  also  the  higher  harmonic  components  of  the  square 
wave  modulated  signal.  The  output  of  this  preamplifier  went 
to  both  the  analog  and  the  digital  acquistion  systems. 

Analog  Data  Acquisition 

The  signal  was  processed  by  the  lock-in  amplifier.  The 
output  went  to  a scaler  (designed  by  Dr.  E.  Voigtman  of  this 
research  lab)  which  acted  as  an  active  divider  for  scaling 
the  10  V output  for  both  the  recorder  and  the  1 V input  of 
the  digital  data  system  which  was  used  for  analyzing  the 
statistics  of  the  analog  output. 

The  reference  signal  was  generated  from  the  reference 
signals  from  the  two  choppers.  A logic  circuit  was  con- 
structed so  that  either  chopper  frequency  could  be  selected, 
or  the  frequency  of  the  sum  or  difference  of  the  two  fre- 
quencies (Appendix  A) . It  was  found  that  the  FFT  of  the 
output  from  the  exclusive  OR  circuit  gave  only  the  sum  and 
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difference  frequencies  when  the  reference  signals  from  both 
choppers  were  combined;  the  frequency  components  of  each  of 
the  single-beam  waveforms  were  suppressed.  This  simple  cir- 
cuit was  followed  by  a high-Q  frequency  selective  amplifier 
which  generated  a bipolar  sinusoidal  reference  signal  from 
the  rectangular  input  signal.  This  frequency  of  this  am- 
plifier was  adjusted  to  be  that  of  the  desired  reference, 
i.e.,  f i , f 2 , or  f^+f ^ . The  second  harmonic  signal  could  be 
easily  measured;  the  lock-in  reference  mode  selector  was 
switched  to  the  2f  mode,  and  the  phase  readjusted  for  maxi- 
mum signal. 

For  the  single  beam  measurements,  the  power  meter  was 
placed  in  the  position  of  chopper  1,  so  that  the  laser  in- 
tensity could  be  continuously  monitored.  The  ouput  of  the 
microvoltmeter  went  to  another  channel  of  the  scaler  . From 
there,  after  scaling,  it  went  to  the  second  channel  of  the 
recorder,  and  could  also  be  sent  to  the  digital  data  an- 
alysis system. 

Digital  Data  Acquisition  and  Signal  processing 

The  computer  system  was  used  for  digital  data  acquisi- 
tion and  processing  using  a fast  Fourier  transform  program 
to  calculate  the  power  spectra  which  gave  the  amplitude  of 
the  various  frequencies  composing  the  signal.  When  using 
the  fast  Fourier  transform  ( F FT)  , there  was  a trade-off 
between  resolution,  frequency  range,  processing  time,  and 
memory  requirements.  The  frequency  of  the  PAR  chopper  could 
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be  varied,  by  changing  the  pulley  ratio  or  the  chopper 
wheel.  The  Rofin  chopper  frequency  could  be  electronically 
varied  continuously  from  several  Hz  to  about  1 kHz.  The 
frequency  of  the  PAR  was  fixed  at  280  Hz,  and  the  Rofin  ar- 
bitrarily set  at  178  Hz.  This  gave  a range  of  frequencies 
of  interest  from  about  100  Hz  ( f 1 — f 2 ) to  about  500  Hz 
( f l + f 2 ^ * The  sampling  frequency  should  be  at  least  twice 
the  highest  frequency  of  interest  (the  Nyquist  frequency)  so 
that  a sampling  frequency  of  at  least  1 kHz  was  required. 

Since  the  signals  were  square  wave  modulated,  the 
higher  harmonics  were  present  and  would  fold  over  if  the 
sampling  frequency  were  only  twice  the  modulation  frequency. 
In  order  to  attenuate  the  higher  harmonics,  a fourth-order 
low-pass  Butterworth  filter  was  built  (Appendix  A)  . With 
the  cut-off  frequency  set  at  700  Hz,  there  was  no  signifi- 
cant attenuation  of  signals  below  500  Hz.  A sampling  fre- 
quency of  4 kHz  was  selected  so  that  the  effect  of  the 
filter  could  be  observed,  and  possible  aliasing  or  fold-over 
minimized.  A 1024-point  transform  gave  2 Hz  resolution 
which  was  adequate  to  resolve  the  various  frequency  com- 
ponents of  the  intermodulated  signal.  The  power  spectrum 
analysis  program  (Appendix  B) , including  the  FFT,  was  based 
upon  one  recommended  by  the  Digital  Signal  Processing  Com- 
mittee of  the  IEEE  Acoustics,  Speech,  and  Signal  Processing 
Society  (25).  For  the  digital  analysis,  8192  data  points 
were  acquired  in  real  time,  at  the  4 kHz  sampling  rate,  and 
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stored  on  disc,  a Hamming  window  of  512  points  was  used  to 
reduce  side  lobes,  and  the  modified  periodograms  averaged. 
A single  1024-point  FFT  took  about  4 s,  or  about  one  minute 
to  perform  one  analysis  (an  average  of  16  FFTs) . The 
resulting  power  spectrum  could  be  printed  out,  or  plotted 
either  using  either  frequency  versus  rms  amplitude/H z1//2  or 
versus  log  power/Hz. 


Experimental  Methods 
Single-Beam  Measurements 

For  initial  characterization  of  the  system  before  going 
to  the  intermodulated  measurements,  a series  of  measurements 
were  made  using  only  the  single-beam  configuration  mentioned 
above.  By  not  using  the  beamsplitter,  approximately  1.2  W 
(over  the  laser  spectral  bandwidth)  could  be  delivered  to 
the  flame  cell  which  should  result  in  a higher  degree  of  sa- 
turation. Two  sets  of  measurements  were  made.  The  first, 
with  constant  power,  but  variable  spot  size,  resulted  from 
adjusting  the  lens  (L2)  so  that  the  beam  was  defocused. 
Thus  the  volume  of  the  flame  illuminated  changed  as  the  lens 
was  moved  from  the  distance  equal  to  the  focal  length. 

In  the  second  set,  the  spot  size  was  kept  nearly  con- 
stant. Here  one  of  the  polarizers  was  placed  in  the  excita- 
tion side  (between  M2  and  M3)  . The  power  meter  was  placed 
at  the  position  of  the  other  chopper  (CHI)  so  that  the 
focused  beam  was  smaller  than  the  aperture  of  the  power 
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meter.  The  laser  was  vertically  polarized  as  a result  of 
the  optics  of  the  dye  laser  cavity  (both  the  dye  stream  and 
the  birefringent  filter  were  at  Brewster's  angle  with 
respect  to  the  laser  beam).  The  power  could  be  varied  con- 
tinuously by  rotating  the  polarizer  from  an  orientation  par- 
allel to  the  incoming  light  to  perpendicular,  which  gave  al- 
most complete  extinction.  The  power  was  varied  from  a maxi- 
mum of  1.2  W to  2 mW  and  the  fluorescence  of  the  1 ppm  Na 
solution  measured.  The  amplitude  of  the  signal  at  the  fun- 
damental and  harmonics  was  determined  from  the  power  spec- 
tra. 


Two-Beam  Measurements 

After  the  single-beam  measurements  had  been  made,  the 
beamsplitter  was  mounted  and  several  series  of  measurements 
performed  to  evaluate  the  two-beam  intermodulation  scheme. 
The  nearly  co-linear  configuration  shown  in  Fig.  3-1  was 
used  so  that  there  would  be  maximal  overlap  between  the  two 
beams.  An  alternative  would  have  been  to  have  the  beams 
crossing  at  an  angle  closer  to  90°  than  180°.  This  would 
have  made  alignment  even  more  critical  and  resulted  in  a 
much  smaller  volume  where  intermodulation  could  occur  . Two 
separate  cases  of  the  co-linear  configuration  were  consi- 
dered: one,  with  the  beams  tightly  focused, and  the  other. 


having  the  beams  defocused  in  the  intersecting  volume. 
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Tightly-focused  beam  alignment 

For  these  measurements,  the  point  of  tightest  focus  was 
chosen.  The  procedure  for  determining  that  point  was  quite 
simple.  One  chopper  frequency  was  selected  as  the  refer- 
ence. The  laser  beam  was  first  coarsely  focused  by  moving 
the  lens  positioner  back  and  forth,  along  the  optical  rail, 
so  that  the  narrowest  portion  of  the  beam  appeared  to  be  in 
the  center  of  the  flame.  The  fluorescence  signal  at  the 
chopper  frequency  was  then  used  as  an  indication  of  the 
alignment.  The  beam  was  moved  vertically  using  the  vertical 
adjustment  of  the  5-axis  lens  positioner.  The  point  of  max- 
imum signal  was  taken  to  be  the  optical  axis.  The  lens 
could  then  be  moved  along  this  axis  using  the  translator  of 
the  lens  positioner  (0.25  in  travel).  For  this  adjustment, 
the  minimum  signal  indicated  the  best  focus,  as  a smaller 
fluorescence  volume  was  illuminated,  and  the  signal  was  sa- 
turated. (Changing  the  vertical  adjustment  would  further 
decrease  the  signal,  but  that  would  be  off  the  optical 
axis.)  Either  the  fundamental  or  the  second  harmonic  signal 
could  be  used  for  this  positioning.  During  this  adjustment 
the  other  beam  need  not  be  blocked.  After  one  beam  was  po- 
sitioned, the  other  beam  would  be  positioned  similarly, 
first  focusing  it  so  it  appeared  to  be  nearly  co-linear  with 
the  first  beam.  After  each  beam  had  been  positioned  so  that 
the  signal  on  the  optical  axis  was  a minimum,  the  reference 
frequency  for  the  lock-in  amplifier  was  changed  from  that  of 
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Measurement  procedures 

For  the  focused  beams,  two  types  of  measurements  were 
made.  In  one,  the  intensity  of  one  beam  was  held  constant, 
and  that  of  the  other  varied  by  placing  a polarizer  between 
the  beamsplitter  and  Ml.  This  beam,  beam  1,  was  the  more 
intense,  due  to  the  higher  reflectivity  of  the  beamsplitter. 
The  polarizer  was  also  used  to  set  the  intensity  of  beam  1 
to  be  the  same  as  that  of  beam  2.  The  total  intensity  was 
varied  by  placing  neutral  density  filters  just  in  front  of 
the  beamsplitter.  The  optical  path  was  not  significantly 
affected,  and  re-alignment  of  the  beams  was  not  necessary  as 
the  neutral  density  filters  were  changed.  This  had  been  a 
problem  when  neutral  density  filters  were  used  instead  of 
the  polarizer  to  vary  the  relative  intensity  of  one  beam. 
For  both  series  of  measurements,  the  fluorescence  signal  was 
measured  using  the  power  spectrum  analysis  program  as  was 
done  for  the  single  beam  measurements. 

For  the  defocused  beams,  the  fluorescence  was  deter- 
mined as  a function  of  the  intensity,  with  the  intensity  of 
the  two  beams  being  nearly  equal,  as  with  the  second  series 
with  the  focused  beams. 

After  the  two  types  of  focusing  had  been  characterized, 
their  use  for  the  reduction  of  scattering  interferences  and 
for  spatial  measurements  was  demonstrated.  Highly  scatter- 
ing matrices  readily  clog  the  type  of  burner  which  was  used 
for  these  studies,  so  an  alternative  was  to  use  the  scatter 
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from  unvaporized  water  droplets  in  the  low  temperature  {^"Ar 
diffusion  flame,  and  also,  as  a worst  case,  the  scatter  from 
water  droplets  with  no  flame.  As  an  example  of  the  spatial 
measurements,  the  flame  was  translated  towards  the  monochro- 
mator, and  away,  to  give  a cross-section  of  the  Na  atom  po- 
pulation across  the  width  of  the  flame.  Both  focused  and 
defocused  beams  were  used.  The  relative  fluorescence  signal 
at  the  sum  frequency  f1+f2  was  recorded  using  the  lock-in 


amplifier . 


CHAPTER  FOUR 
RESULTS  AND  DISCUSSION 


All  power  spectra  plotted  are  averages  of  16  1024  point 
transforms  as  indicated  in  Chapter  3 and  Appendix  C,  with  a 
resolution  of  2 Hz.  This  corresponds  to  a two-second  time 
series.  The  ordinate  for  the  linear  plots  is  scaled  propor- 
tional to  rms  amplitude/Hz1/2  while  that  for  the  semi-log 
plots  is  the  power  spectral  density,  amplitude2/Hz , in  deci- 
bels. The  linear  plots  are  useful  for  qualitative  compari- 
son of  signal  levels,  i.e.,  signals  of  similar  order  of  mag- 
nitude. The  semi-log  plots  reveal  more  of  the  fine  detail 
of  the  spectral  composition  of  the  modulated  signal  and 
gives  a better  understanding  of  the  interactions. 


Single-Beam  Modulation 

Typical  power  spectra  for  square-wave  modulation  are 
given  in  Figures  4-1  and  4-2.  The  magnitude  of  the  funda- 
mental, the  second  and  the  third  harmonic  are  plotted  versus 
the  measured  power  in  Figure  4-3.  The  signals  plotted  are 
taken  from  the  print-out  of  the  power  spectra  and  correspond 
to  the  peak  value  with  the  baseline  subtracted.  Similar 
results  are  obtained  by  integrating  the  data  corresponding 
to  the  given  peak. 
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Figure  4-1.  Typical  flame  fluorescence  power  spectrum  for 
1 ppm  Na,  with  single  beam  focused  at  high  power  (425  mw) . 
Chopper  frequency  280  Hz. 

a)  Linear  plot 

b)  Semi-log  plot 
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FREQUENCE  (llz.) 


Figure  4-2.  Typical  power  spectrum  for  fluorescence, 

1 ppm  Na,  single  beam,  square-wave  modulated,  focused 
at  lowpower  (3.5  mW) 

a)  Linear  plot 

b)  Semi-log  plot.  Second  harmonic  is  significantly 
reduced  compared  to  Fig  4.1 


Amplitude  /Hz 
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a) 


b) 


Frequency  (Hz) 


Figure  4-3.  Fluorescence  signal  at  the  fundamental,  second 
and  third  harmonics  vs  laser  power.  Single  focused 
beam,  square-wave  modulated. 

• Fundamental 
■ Second  harmonic 
A Third  harmonic 


rms  amplitude  Hz-I/2  (a.u.) 
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The  magnitude  of  the  second  harmonic  was  unexpected. 
Since  the  square  of  a square  wave  is  another  square  wave,  it 
might  be  expected  that  the  two  would  have  the  same  frequency 
content.  Freuholz  and  Gelbwachs  (13)  found  a fundamental  to 
second-harmonic  ratio  of  about  6 for  square-wave  modulation 
with  intensities  between  3 and  7 kW/cm2.  Here  the  ratio  is 
about  3.  From  the  first-order  approximation  for  square-wave 
modulation,  a factor  similar  to  Freuholz' s would  be  expect- 
ed. By  examining  the  Fourier  series  expansion  of  a square 
wave,  Eqn.  2.13,  and  the  square  thereof, 

2 4 1 i 9 

f(fc)  = ~2  (coscot COS3  cot  + -COS5  tot  ...) 

77  3 5 

one  can  see  the  origin  of  a second-harmonic  component,  from 

9 

the  identity  cos  A = (l+cos2A)/2. 

The  shape  of  the  curve  is  typical  of  saturation  curves 
where  the  laser  beam  is  focused  to  a spot  that  is  small  com- 
pared to  the  flame  thickness.  Only  where  the  spatial  dis- 
tribution of  the  laser  irradiance  is  uniform  in  the  observed 
part  of  the  flame  is  a plateau  found  as  in  Fig.  2-lb  (24) 

It  should  be  noted  that  the  point  corresponding  to  the 
highest  power,  520  mw  average  (with  chopper  in  place)  or 
1.04  W peak,  was  the  only  one  with  no  polarizer  in  the  beam 
path,  so  the  spatial  distribution  might  have  been  somewhat 
d if ferent . 
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Single-beam  modulation  has  the  advantages  of  requiring 
fewer  optical  components  and  having  less  critical 
adjustments.  Spatial  resolution  is  determined  by  the  volume 
of  the  atom  cell  that  is  saturated  plus  the  projection  of 
that  upon  the  monochromator  entrance  slit.  Thus  optical 
alignment  may  be  critical  if  spatial  information  is  desired. 
The  magnitude  of  the  second  harmonic  will  be  considered  in 
the  evaluation  of  the  two-beam  system,  to  compare  the  effec- 
tiveness of  harmonic  mixing  versus  intermodulation. 

Double-Beam  Modulation 
Tightly-Focused  Beams 

The  complexity  of  the  intermodulation  power  spectrum  is 
shown  in  Figure  4-4  where  the  major  intermodulation  products 
are  labeled.  Figure  4-5  gives  two  typical  power  spectra  for 
the  double-beam  measurement.  in  this  first  series  of  meas- 
urements, the  intensity  of  beam  2 (passing  through  chopper 
2)  was  held  constant  and  the  intensity  of  beam  1 was  attenu- 
ated, so  that  it  varied  from  almost  twice  the  other  to  al- 
most complete  extinction.  The  peak  intensity  of  beam  1 was 
nominally  320  mW.  The  results  are  shown  in  Fig.  4-6.  The 
signal  at  f1+f2  is  maximal  closer  to  240  mW  which  may  be  at- 
tributed to  imperfect  matching  of  the  two  beams;  the  ratio 
of  the  fundamental  signals  (^/f^  were  1.15  and  1.08  with 
320  m'W  and  240  mW  excitation,  indicating  there  was  a closer 
match  in  intersecting  intensities  at  the  240  mW  level.  This 
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4-5.  Power  spectra  of  intermodulated 
fluorescence  signal  for  tightly-focused 
Beam  2,  320  mw 

a)  Beam  1,  high  intensity 

b)  Beam  2,  low  intensity 
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a) 


Frequency  (Hz) 


Figure  4-6.  Comparison  of  second-harmonic  and  sum 
frequency  signals  with  square-wave  modulation 
and  tightly-focused  beams.  Intensity  of  one  beam 
held  constant.  (320  mw) 


rms  amplitude  Hz-1'*  (a.u.) 
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curve  indicates  there  is  an  optimum  for  the  relative  inten- 
sities of  the  two  beams,  as  expected  from  theory 
(Eqn.  2.18),  and  that  experimentally,  approximately  equal 
intensities  are  required  for  the  highest  intermodulation 
s ignal . 

Since  the  signal  at  2*f^  is  expected  to  be  a function 
of  the  square  of  I1  while  the  signal  at  (f1  + f2)  is  a func- 
tion of  I^*I2»  with  being  held  constant,  it  is  not 
surprising  that  the  second  harmonic  component  becomes  larger 
than  the  intermodulation  signal  as  ^ becomes  larger  than 
I2*  There  appears  to  be  less  than  a factor  of  two  differ- 
ence between  the  second  harmonic  and  the  intermodulation 
signal.  The  single-beam  signal  might  very  well  be  larger 
than  that  from  the  dual-beam  system,  as  the  laser  beam  is 
not  split,  nor  suffer  the  losses  from  the  additional  optical 
components.  Thus,  square-wave  modulated  harmonic  saturation 
spectroscopy  should  be  comparable  with  the  double-beam  in- 
termodulated  spectroscopy.  it,  however,  has  limited  poten- 
tial for  spatial  resolution. 

Defocused  Beams 

The  differences  between  the  intermodulation  and  harmon- 
ic products  for  the  focused  and  defocused  beams  are  shown  in 
Figs.  4-7  and  4-8.  The  tightly-focused,  intersecting  beams 
have  a high  second-harmonic  content,  which  is  to  be  expected 
from  the  single-beam  results.  Since  the  volume  viewed  by 
the  monochromator  is  much  larger  than  the  intersection 


Figure  4-7.  Comparison  of  intermodulation  signals 
from  beams  tightly  focused  at  intersection  volume 
and  defocused  at  intersection  volume.  Laser 
power  360  mw.  Linear  power  spectra. 

a)  Tightly-focused  beams,  gain  of  200 

b)  Defocused  beams,  gain  of  20 
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Frequency  (Hz) 


Figure  4-8.  Comparison  of  intermodulation  signals 
from  beams  tightly  focused  at  intersection  volume 
and  defocused  at  intersection  volume.  Laser 
power  360  mw.  Semi-log  plot  of  power  spectra. 

a)  Tightly-focused  beams,  gain  of  200 

b)  Defocused  beams,  gain  of  20 
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volume  of  the  two  beams  in  the  case  of  the  tightly-focused 
beams.  Fig.  4-7a,  it  is  not  surprising  that,  in  this  case, 
the  single-beam  components,  e.g.,  signals  at  and  2f±f  are 
large  compared  to  the  intermodulation  signal  at  f^f  For 
the  defocused  case,  the  signal  levels  of  the  fundamental  are 
generally  5-10  times  larger  than  those  with  the 
tightly-focused  beams  at  the  same  power  levels,  due  in  part 
to  the  larger  illuminated  volume,  with  the  defocused  beams, 
the  second  harmonic  signal  is  not  discernible,  and  the  in- 
terraodulation  is  relatively  large  compared  to  the  focused 
double-beam  intermodulation  signals. 

The  effect  of  varying  the  laser  power  for  the  defocused 
beams  is  shown  by  Fig.  4-9,  where  the  power  of  the  two 
beams  is  nominally  the  same.  The  fundamental  frequency  com- 
ponent is  leveling  off  with  increasing  power  while  the  in- 
termodulation product  continues  to  increase.  Thus,  for  high 
resolution,  one  could  use  the  two  focused  beams,  while  for 
good  scattering  rejection,  one  would  get  larger  signals  and 
average  over  a greater  volume  using  the  defocused  beams. 

Scattering  Rejection 

The  effectiveness  of  the  intermodulation  in  reducing 
scatter  is  shown  in  Figs.  4-10  and  4-11,  where  the  power 
spectrum  for  scatter  are  plotted,  for  scatter  from  water 
droplets  in  a ^"Ar  diffusion  flame  (no  C^)  and  from  water 
droplets  with  no  flame.  in  neither  case  is  there  any  com- 
ponent at  the  sum  or  difference  frequencies  f^+/-f 


The 


Figure  4-9.  Fluorescence  signal  at  fundamental  and 
sum  frequency  for  two-beam  square-wave  modulation. 
Beams  nearly  equal  intensity  and  defocused  at 
point  of  intersection;  total  power  varied. 


rms  amplitude  Hz-,/2  (a.u.) 
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aspir ated 
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10.  Power  spectrum  for  scatter  from  water 
in  H--Ar  diffusion  flame.  Gain  = 20,  two 
beams,  360  mw  each, 
linear  plot 
semi-log  plot 
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a) 


Frequency  (Hz) 


b) 


Frequency  (Hz) 


Figure  4-11.  Power  spectrum  for  scatter  from  water 
aspirated  with  flame  off.  Gain  = 10  , two  defocused  beams, 
360  mw  each 

a)  linear  plot 

b)  semi-log  plot 
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shot  noise  from  the  large  scatter  signal  gives  the  nearly 
uniform  density  out  to  about  700  Hz,  the  cut-off  frequency 
of  the  Butterworth  filter.  The  roll-off  of  this  white  noise 
power  spectrum  indicates  that  signals  above  1.5  kHz  are  at- 
tenuated by  a factor  of  at  least  10-3  (60  db) . 

The  reduction  in  the  scattering  signal  may  be  shown 
more  clearly  in  the  time  domain,  comparing  lock-in  amplifier 
output  signals  at  the  fundamental  frequency  and  at  the  in- 
termodulation frequency,  as  in  Fig.  4-12.  This  represents 
a significant  reduction  in  the  scatter  signal  and  indicates 
the  potential  of  using  intermodulation  saturation  spectros- 
copy for  the  reduction  of  scatter  interferences. 

Spatial  Resolution 

As  an  example  of  the  use  of  the  intermodulation  signal 
for  spatial  resolution,  the  Na  atom  density  in  the  flame  was 
mapped  by  translating  the  flame  towards  the  monochromator 
and  back  to  get  the  distribution  across  the  flame  width. 
The  flame  was  moved  in  approximately  1 mm  steps.  The  length 
of  the  volume  element  measured  was  determined  by  either  the 
volume  in  which  intermodulation  saturation  was  occurring,  or 
by  the  2 mm  slit  height  (slit  oriented  parallel  to  the 
beam),  whichever  was  smaller.  The  profile  was  measured  with 
the  focused  beams,  and  also  the  defocused  beams.  The  defo- 
cused  beams  gave  approximately  10  times  larger  a signal,  and 
had  better  signal/noise  ratio,  but  either  was  adequate  for 
this  purpose.  The  scan  across  the  flame  shown  in  Fig.  4-13 


Figure  4-12.  Recorder  tracing  of  scatter  signal, 
defocused  square-wave  modulated  beams,  360  mw  each 

a)  Hydrogen  diffusion  flame,  as  in  Fig.  4-10, 
lock-in  amplifier  reference  frequency  f, 

b)  as  in  a,  f ^ reference  frequency  1 

c)  argon  + water  aroplets , as  in  4-11.  Signal 
measured  at  f,  . 

d)  as  in  c,  signal  measured  at  f^  + f2 
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Figure  4-13.  Sodium  atom  distribution  across  width  of  flame 
1 ppm  Na  in  H^-O^-Ar  flame. 

a)  Schematic  diagram;  burner  moved  along  optical  axis 
of  monochromator,  ~ 1 mm  between  observations 

b)  intermodulated  fluorescence  signal  measured  with 
lock-in  amplifier  reference  frequency  of  f. +f ~ , 
using  defocused,  nearly  co-linear,  intersecting 
beams 
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gives  the  relative  Na  atom  distribution  2 cm  above  the 
of  the  burner  head,  under  the  flame  conditions  used  in 
s tudy . 


top 
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CHAPTER  FIVE 
CONCLUSIONS 


The  intermodulation  component  of  the  fluorescence  sig- 
nal in  atomic  fluorescence  may  be  a significant  fraction  of 
the  saturation  signal.  By  using  two  beams  modulated  at  dif- 
ferent frequencies,  a resultant  sum  or  difference  frequency 
may  be  selected  to  avoid  frequency  regions  with  high  noise 
components.  Harmonic  saturation  spectroscopy  is  a special 
case  of  intermodulated  saturation  spectroscopy  in  which  the 
frequency  and  position  of  the  two  beams  are  the  same.  The 
same  theoretical  considerations  apply  in  determining  the  mo- 
dulation products. 

The  use  of  intermodulated  atomic  fluorescence  for  the 
reduction  of  scattering  and  for  spatial  resolution  has  been 
demonstrated  using  Na  in  a H2-02-Ar  flame.  if  saturation 
can  be  achieved,  the  signal/noise  ratio  may  be  improved  by 
defocusing  the  beams  and  illuminating  a larger  volume.  If 
fine  spatial  resolution  is  required,  the  beams  may  be  tight- 
ly focused  using  short  focal  length  lenses.  Sodium  has  been 
used  extensively  as  a seed  for  studies  of  flame  chemistry. 
It  is  of  interest  for  flame  diagnostics,  combustion  diagnos- 
tics, and  as  a model  from  which  techniques  have  been  de- 
veloped for  the  study  of  radicals  in  combustion  processes, 
as  well  as  for  physical  studies. 
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Farther  application  of  intermodulated  flame  fluores- 
cence would  require  the  use  of  lasers  which  can  saturate 
other  elements  of  analytical  interest.  Pulsed  dye  lasers 
have  much  higher  peak  powers  and  cover  a wider  range  of 
wavelengths.  For  such  devices,  a suitable 
modulation/detection  scheme  would  need  to  be  developed, 
probably  somewhat  more  elaborate  than  the  preamplifier  and 
phase-sensitive  amplifier  or  analog- to-d ig i tal  converter 
used  here  for  the  modulated  cw  dye  laser. 


APPENDIX  A 

CIRCUITS  FOR  INTERMODULATION  MEASUREMENTS 
Lock-in  Amplifier  Reference  Circuit 
The  circuit  shown  in  Fig.  A-l  was  constructed  to  gen- 
erate a signal  that  could  be  used  as  the  reference  waveform 
for  the  intermodulation  signal.  input  and  output  waveforms 
are  shown  in  Fig.  A-2.  The  power  spectra  of  the  ouptputs 
are  shown  in  Figs.  A-3  and  4.  The  exclusive-OR  circuit 
generates  only  the  sum  and  difference  frequencies  (and  other 
higher  frequency  components  of  those  waveforms).  Thus,  the 
requirements  for  a frequency  selective  amplifier  to  give  a 
single  frequency  output  are  not  so  demanding,  using  this 
configuration.  Schmitt  trigger  gates  were  used  as  the 
chopper  reference  signals  are  not  necessarily  TTL  compati- 
ble. The  operational  amplifiers  are  used  in  a non-inverting 
voltage-follower  configuration  as  output  buffers  for  analog 
signals.  A frequency-selective  amplifier  (par  210)  was  used 
to  select  a particular  frequency  component  of  the  output 
waveform,  and  to  generate  a bipolar  sine  wave  which  was  used 
as  the  reference  signal  for  the  lock-in. 

The  alternate  output,  from  the  NAND  gates,  was  used  for 
test  purposes  as  it  simulates  the  fluorescence  signal  in  the 
case  of  complete  saturation. 
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Figure  A-2.  Lock-in 

a)  Reference 

b)  Reference 

c)  NAND  gate 

d)  Exclusive 


amplifier  reference 

1 waveform,  f, 

2 waveform,  f^ 
output,  B 

OR  output,  A 


circuit  waveforms. 
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Figure  A-3.  Power  spectrum  of  exclusive  OR  reference 
output  waveform 

a)  Linear  plot 

b)  Semi-log  plot 


s Amplitude  / Hz 
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a) 


b) 


Frequency  (Iiz) 


Figure  A-4.  Power  spectrum  of  NAND  gate  reference 
output  waveform 

a)  Linear  plot 

b)  Semi-log  plot 


s Amplitude  / H2 
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Anti-aliasing  Filter 

A low-pass  filter  was  constructed  to  prevent  aliasing, 
i.e.,  higher  frequency  signals  appearing  at  lower  frequen- 
cies, with  the  FFT.  In  order  to  obtain  a sharp  cutoff,  a 
fourth-order  Butterworth  filter  was  constructed  using  a di- 
gital switched-capac ito r filter  (National  MF10).  The 
Butterworth  filter  has  maximally  flat  amplitude  response  in 
the  passband,  though  it  has  a non-linear  phase  characteris- 
tic which  results  in  over-shoot  in  the  time  domain  (19),  but 
this  was  not  considered  significant  for  the  measurement  of 
power  spectra. 

The  ratios  R2A/R2B  and  R3A/R3B,  which  determine  the 
filter  characteristics,  were  obtained  by  trimming  available 
components  to  give  the  desired  frequency  response.  The 
resistors  R01  and  R02  were  selected  to  give  minimum  DC 
offset  of  the  internal  amplifiers,  and  thus  maintain  maximum 
dynamic  range.  The  overall  voltage  gain  was  approximately 
one  for  signals  within  the  pass-band.  The  maximum  input 
voltage  should  be  less  than  5.6  Vp-p  for  this  design,  ac- 
cording to  the  manufacturer's  application  notes.  Since  the 
A/D  converter  with  which  this  filter  was  used  had  a range  of 
+/-  1 V,  the  input  signal  was  kept  less  than  2 Vp-p. 

This  switched-capacitor  filter  building  block  uses  an 
external  clock  to  set  the  cutoff  frequency,  which  makes  it 
quite  versatile  for  use  with  the  FFT  for  power  spectrum  an- 
alysis. Analog  filters,  especially  higher-order  filters, 
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require  several  precisely  matched  components  for  each  choice 
of  frequencies.  This  filter  requires  precise  matching  only 
of  the  resistors  to  estaolish  the  filter  characteristics, 
and  then  a signal  generator  can  be  used  to  vary  the  cutoff 
frequency . 

The  amplitude  response  was  checked  using  a signal  gen- 
erator as  an  input  and  determining  the  amount  of  attenua- 
tion. in  order  to  have  minimal  attenuation  over  the  fre- 
quency range  of  interest,  a cutoff  frequency  of  700  Hz  was 
selected.  This  gave  a 2%  attenuation  at  500  Hz.  The  active 
low-pass  filter  at  the  output  reduces  some  of  the  clock  and 
switching  noise  which  is  characteristic  of  this  type  of  dev- 
ice. it  also  acts  as  a buffer  amplifier  to  avoid  any  load- 
ing effects  on  the  filter  output. 


APPENDIX  B 

ESTIMATION  OF  POWER  SPECTRA 

The  frequency  components  of  the  fluorescence  signal 
were  determined  by  estimating  the  power  spectra  using  the 
FFT.  The  program  was  based  on  one  recommmended  by  a Digital 
Signal  Processing  Committee  of  the  IEEE  Acoustics,  Speech 
and  Signal  Processing  Society  (25).  They  had  tested  the 
program  on  various  computers,  including  a DEC  POP  11/34A 
similar  to  the  one  used  here,  and  found  it  suitable  to  re- 
commend for  use.  The  modifications  made  here  did  not  in- 
volve changing  the  algorithm,  but  only  the  input  and  output 
for  use  with  this  real-time  data  acquisition  system. 

The  program  is  based  on  the  modified  periodogram  method 
of  Welch  (26),  which  is  computationally  efficient  and  gives 
good  results  (27).  The  periodogram  is  a measure  of  the  rms 
amplitude  of  a short  data  record.  The  periodogram  is  modi- 
fied through  the  use  of  a window,  which  compensates  for  the 
finite  data  record  and  reduces  side-lobes  on  the  resultant 
periodogram  (28),  as  shown  in  Figs.  b-1  and  B-2 . A number 
of  modified  periodograms  are  averaged  to  give  an  estimate  of 
the  power  spectrum.  Several  other  FFT  subroutines  from  the 
same  reference  were  evaluated  to  speed  the  data  analysis  (4 
seconds  for  a 1024  point  complex  transform) . no  significant 
improvement  (less  than  a factor  of  two,  accompanied  by  heavy 
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Figure  B-l.  Power  spectrum  of  sine  wave  input. 
Rectangular  window,  otherwise  same  conditions  as 
experimental  measurements 

a)  linear  plot 

b)  semi-log  plot 
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Figure  B-2.  Power  spectrum  of  sine  wave  input. 
Hamming  window;  same  conditions  as  experimental 
measurements . 

a)  linear  plot 

b)  semi-log  plot 
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utilization  of  storage)  was  obtained  with  floating  point 
routines.  The  integer  routine  in  the  DEC  Lab  Applications 
Package  was  a factor  of  four  faster  but  was  not  compatible 
with  the  use  of  a window  for  the  data. 

The  A/D  subroutine  (originally  written  by  d.L.  Bolton 
and  modified  for  use  under  the  TSX  operating  system)  could 
acquire  data  at  a rate  up  to  44  kHz.  The  maximum  length  re- 
cord was  8192  (8  k)  data  points,  due  to  the  size  of  the  For- 
tran i/o  package,  a number  of  8 k records  could  be  stored 
on  disc  and  subsequently  analyzed  with  the  power  spectrum 
analysis  program.  Results  could  be  plotted  on  a video 
graphics  terminal  (DEC  VT100  with  Matrox  GT-600  graphics)  or 
with  a digital  plotter  (Houston  Instruments  DMP-4)  . a 
graphics  subroutine  package  from  the  DEC  Users  Group 
(SFGL70)  was  used  for  the  video  plotting.  a set  of  su- 
broutines (HIPLIB)  written  by  D.L.  Bolton  (29)  was  used  for 
the  hard-copy  output. 
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